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Abstract: The beam-helicity asymmetry in exclusive electroproduction of real photons by the 
longitudinally polarized Hera positron beam scattering off an unpolarized hydrogen target is mea- 
sured at Hermes. The asymmetry arises from deeply virtual Compton scattering and its interfer- 
ence with the Bethe-Heitler process. Azimuthal amplitudes of the beam-helicity asymmetry are 
extracted from a data sample consisting of ep — > ep"f events with detection of all particles in the 
final state including the recoiling proton. The installation of a recoil detector, while reducing the 
acceptance of the experiment, allows the elimination of background from ep —J- eNn'y events, which 
was estimated to contribute an average of about 12% to the signal in previous Hermes publications. 
The removal of this background from the present data sample is shown to increase the magnitude 
of the leading asymmetry amplitude by 0.054 ± 0.016 to -0.328 ± 0.027 (stat.) ± 0.045 (syst.). 
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1 Introduction 

Generalized Parton Distributions (GPDs) [1-3] describe the soft, non-pcrturbativc part of hard 
exclusive reactions, e.g., hard exclusive leptoproduction of a real photon or a meson leaving the 
nucleoli intact (possibly modulo isospin rotation). These distributions have quickly increased in 
importance in QCD spin physics since it was shown that they can provide access to the total angular 
momentum carried by quarks in the nucleon (and also by gluons, in principle) [4]. The resulting 
intense theoretical activity is exemplified by the demonstration that GPDs can be considered as form 
factors, dissected in longitudinal nucleon momentum, describing transverse density distributions of 
quarks (and gluons) [5] ("nucleon tomography"). Hard exclusive processes provide experimental 
access to GPDs. One of the experimental challenges in these measurements is the selection of truly 
exclusive final states, discriminating against the excitation of baryonic resonances. One solution to 
this problem is the detection of all particles in the final state with adequate kinematic resolution. 

Generalized parton distributions depend on four kinematic variables: t, x, £, and Q 2 . The 
Mandelstam variable t = (p — p') 2 is the square of the difference between the initial (p) and final 
(p') four-momenta of the target proton. The variable x is the average of the initial and final 
fractions of the (large) target longitudinal momentum that is carried by the struck parton, and the 
variable £, known as the skewness, is half of the difference between these fractions. The evolution 
of GPDs with the photon virtuality Q 2 = — q 2 is analogous to that of parton distribution functions, 
with q = k — k' being the difference between the four-momenta of the incident and the scattered 
leptons. Currently, there exist no hard exclusive measurements that provide access to x. Because 
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of the lack of consensus about how to define £ in terms of experimental obscrvables, the results are 
typically reported by Hermes as projections in xb = Q 2 /(2pq), to which £ can be related through 
£ ~ xb/(2 — xb) in the generalized Bjorken limit of large Q 2 and fixed xb and t. 

Several GPDs describe various possible helicity transitions of the struck quark and/or the 
nucleon as a whole. At leading twist (i.e., twist-2) and for a spin-1/2 target such as the proton, four 
chiral-even GPDs (H, H, E, E) are required to describe processes that conserve the helicity of the 
struck quark. Deeply Virtual Compton Scattering (DVCS), i.e., the hard exclusive leptoproduction 
of a real photon, has the most reliable interpretation in terms of GPDs among all presently practical 
hard exclusive probes. The measurement of the DVCS process on unpolarized protons is most 
sensitive to GPD H, which describes the transition that conserves the helicities of both the struck 
quark and the nucleon. Such measurements were performed by Hermes [6-9], HI [10-13], and 
Zeus [14, 15] at Hera, by the Hall A Collaboration [16] and by Clas [17-19] at Jefferson Lab. 

This paper reports a kincmatically complete measurement of DVCS for a polarized lcpton 
beam on an unpolarized hydrogen target with detection of all particles in the final state, including 
the recoil proton. It is the first measurement of this kind reported by Hermes. The results of 
this measurement are compared to measurements without detection of the recoil proton, while 
accounting for the difference in acceptance in the two techniques. Measurements without recoil- 
proton detection have also been reported in previous publications. 



2 Constraining GPDs through DVCS 

The four- fold differential cross section for exclusive single-photon production, ep —> ep-f, on an 
unpolarized proton target is given by [20] 



d 4 <7 xbc 6 
&Q 2 dxBdt&4> = 32(2^)4Q4 V / 1 = 



\Tep^fepj\ 7 (2-1) 



where Tep^ep^ is the scattering amplitude for this process. Here, e is the elementary charge and 
e = 2xb^q- with M p the proton mass. The angle (j> denotes the azimuthal orientation of the photon 
production plane with respect to the lepton scattering plane, as indicated in figure 1. This definition 
follows the Trento conventions [21]. 

In DVCS, a real photon is radiated from the struck parton, see figure 2(a). There is an- 
other process that contributes to the channel ep — > epj: the Bethe-Heitler (BH) process, where a 
bremsstrahlung photon is radiated from the incident or scattered lcpton, sec figure 2(b). The DVCS 
and BH processes have the same initial and final states, and therefore their scattering amplitudes 
interfere: 

\%ij~ie P f\ 2 = |7bh| 2 + |7dvcs| 2 +Z, (2-2) 

with the term X representing the interference between the scattering amplitudes of the BH process, 
Tbh, and the DVCS process, 7dvcs : 

I = TbbTbvcs + 7dvcs7bh- (2-3) 

The BH scattering amplitude is calculable in QED using the form factors of the proton measured 
in elastic scattering. 

The contributions to the cross section expressed through eqs. (2.1) and (2.2) can each be 
expanded in a harmonic series with respect to the angle <f> [20]: 
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Figure 1. Momenta and azimuthal angle for exclusive electro-production of real photons in the target 
rest frame. The quantity <f> denotes the angle between the lepton scattering plane containing the three- 
momenta k and k' of the incoming and outgoing lepton, respectively, and the plane correspondingly defined 
by q — k — k' and the three-momentum q' of the real photon. Also indicated are the polar angle (9 7 * 7 
between the three-momenta of the real (q') and virtual (q) photons, and the three-momentum of the recoil 
proton (p 1 ). 




irB » |2= ^^»( c » H+ S c » Hcos H' (2 - 4> 

itdvci 2 = jw s (cr cs +Ecr as » s w)+A S r cs ™.») . (2-5) 

n=l / 
3 2 \ 

Co + Yl c » C0S M) + A s n sin M) • (2-6) 



-eiK x 



Here, Vi((f>) and Vq.{<\>) are the lepton propagators for the BH process, and A = ±1 and ei = ±1 are 
respectively the helicity and unit charge of the beam lepton. The quantities -£Tbh = 1/( x b^(1 + £2 )) 2 > 
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-K'dvcs = 1/Q 2 7 an d Kx = l/(ytxB) are kinematic factors independent of 0, with y = (pq)/(pk). 
All beam-helicity dependent terms enter with sinusoidal harmonics due to parity conservation. 

At the kinematic conditions of Hermes, the square of the DVCS scattering amplitude yields 
only a small contribution to the ep — > ep7 cross section, while the square of the BH scattering 
amplitude is much larger. Therefore, the contribution of the DVCS scattering amplitude to the 
cross section enters mainly through its interference with the BH scattering amplitude, thereby 
giving rise to cross-section asymmetries with respect to beam charge and beam helicity (and/or 
target polarization). For a longitudinally polarized beam and an unpolarized target, the cross 
section of eq. (2.1) can be expressed as: 



d a 



dQ 2 d:EB dti 



clu(0, e e , A) = cruu(0, e-i) [1 + A.A L u(0, e-e)} 



(2.7) 



where ctlu (fuu) denotes the differential cross section for longitudinally polarized (unpolarized) 
beam and unpolarized target. On the right-hand side, the other kinematic dependences are omitted 
for brevity. For an unpolarized beam, the cross section reads: 



cuu(0, ee) 



x B ef 



32(27r) 4 Q 4 v / IT? 



'^DVCS C 



DVCS 



^ BH (cr+ELi^ H cos( 



E 



c° vcs cos(n 




e e K x (cj + ELi4 C0S ( 



Vi{<t>)v 2 {4>) 

In eq. (2.7), A\ J \j{4> 1 ei) denotes the single-charge beam-helicity asymmetry: 



A 



IAM 



gLu(^, e e , A = +1) - QLuQft, e e , A = -1) 
, A = +1) + cr L u(0, ee, A = -1) 



Kovcs s? VCS sin< 



: sin(n0) 



(2.9) 
(2.10) 



It can be seen that Aiaj(4>, ei) is a mixture of contributions from beam-charge dependent interfer- 
ence and beam-charge independent squared DVCS terms in both its numerator and denominator, 
the latter being identified with o~uu- For data collected with only one beam charge, the two terms 
in eq. (2.10) containing sf v and sf cannot be disentangled. However, at Hermes kinematic con- 
ditions, the asymmetry is expected to be dominated by the term containing s\ [20]. Measurements 
that disentangle the contributions to the cross section from the interference term and the squared 
DVCS term of eq. (2.2) confirm this expectation by finding that the asymmetry related to the latter 
(and hence sf vcs ) is negligible compared to that arising from the former [8]. 



The leading-twist Fourier coefficient sj is related to GPDs via a complex function C„ 



y(2 - y) 3m(C. 



1 ) 

unp / 



(2.11) 



(The coefficient sf enters at higher twist.) The function Cj np is a linear combination of the Compton 
Form Factors (CFFs) H, H, and £ [20], which are flavor sums of convolutions of the corresponding 
GPDs H , H, and E with hard scattering coefficient functions. This linear combination reads: 



^unp — F\H 



2-xb 



(Fi+F 2 )n-—F 2 £ 



4A/2- 



(2.12) 
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where Fi and F2 are respectively the Dirac and Pauli form factors of the nuclcon. At Hermes 
kinematic conditions, where both xb and —t/Mp arc of order 0.1, to first approximation the con- 
tributions from CFFs H and £ are negligible in eq. (2.12) with respect to that of T-L since they are 
kincmatically suppressed by an order of magnitude or more. Therefore, in this approximation, the 
behavior of Cjf np is determined by the CFF H, and hence the GPD H . Thus, GPD H is constrained 
by measurements of the beam-hclicity asymmetry Auj. 

3 The Hermes experiment in 2006-2007 

During the Hera winter shutdown 2005/2006, a recoil detector [22] was installed in the target 
region of Hermes. The configuration of the forward spectrometer [23] remained unchanged. The 
main purpose of the recoil detector has been the detection of the recoil target proton in order to 
enhance access to hard exclusive processes at Hermes, in particular to DVCS. Data were collected 
from the recoil detector in conjunction with the forward spectrometer in 2006 and 2007 using the 
Hera electron or positron beam of energy 27.6 GeV scattering off a target of unpolarizcd hydrogen 
or deuterium gas internal to the Hera lepton storage ring at Desy. 

The Hera lepton beam was transversely self-polarized by the emission of synchrotron radiation 
[24]. Longitudinal polarization of the beam in the target region was achieved by a pair of spin 
rotators located upstream and downstream of the experiment [25] . The sign of the beam polarization 
was reversed three times over the running period. Two Compton backscattering polarimeters [26, 27] 
independently measured the longitudinal and transverse beam polarizations. 

The commissioning of the recoil detector was completed in 2006 after the switch of the Hera 
lepton beam from electrons to positrons. Therefore, for the analysis of the beam-hclicity asymmetry 
considered here, data collected with only one lepton beam charge but both beam-helicity states are 
available. For this data set, the average beam polarization was Pi — 0.402 (—0.394) for positive 
(negative) beam helicity, with the total relative uncertainty of 1.96% [28]. 

The scattered lepton and particles produced in the polar-angle range 0.04 rad < 9 < 0.22 rad 
were detected by the Hermes forward spectrometer. The average lepton-identification efficiency 
was at least 98% with hadron contamination of less than 1%. The produced particles emerging 
at large polar angles and with small momenta were detected by the Hermes recoil detector in 
the polar-angle range 0.25 rad < 8 < 1.45 rad, with an azimuthal coverage of about 75%. The 
lower- momentum detection threshold for protons was 125 MeV for this analysis. 

The recoil detector surrounded the Hermes target cell and consisted of several subcomponents 
embedded in a solenoidal magnetic field with field strength of 1 T. A detailed description of the 
recoil-detector components is given in ref. [22]. Figure 3 shows a schematic view. 

The innermost active detector component, surrounding the target cell, was a Silicon Strip 
Detector (SSD) made up of two layers. Each layer consisted of eight double-sided sensors with 
orthogonal strips on opposite sides. In order to minimize the momentum threshold for proton 
detection, the amount of passive material between the 75 y,m thick target-cell wall and the first 
layer of sensors was minimized by placing the SSD and the front-end read-out electronics inside the 
Hera beam vacuum as close as 5 cm from the electron or positron beam. Each sensor had an area 
of 9.9 cm x 9.9 cm and a strip pitch of 758.2 fim with individual sensor thicknesses varying between 
295 /im and 315 /im. 

Protons with momenta larger than 250 MeV passed through the 1mm thick wall of the alu- 
minum vacuum chamber and were detected by the Scintillating Fiber Tracker (SFT). The latter 
consisted of two coaxial barrels of scintillating fibers of 1 mm diameter, the inner barrel at a radius 
of about 11.5 cm and the outer barrel at a radius of about 18.5 cm. The active length of both barrels 
was 28 cm. A barrel was made of two adjacent sub-barrels each consisting of two layers of fibers. 
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PD 




Figure 3. Schematic diagram of the Hermes recoil detector (cross-section view). The Hera lepton beam 
is perpendicular to the paper plane. The cross section of the target cell is shown as ellipse. The tracking 
layers are indicated, from inside to outside: inner and outer layers (in diamond shape) of the Silicon Strip 
Detector (SSD), inner and outer barrels (circles) of the Scintillating Fiber Tracker (SFT). Space-points are 
indicated by crosses. The SSD modules are located inside the vacuum chamber (dashed circle). The Photon 
Detector (PD) shown as a dash-dotted circle is not used in the present analysis. The magnet (not shown) 
surrounds the detector assembly. Also shown are examples of tracks reconstructed from two, three, and 
four space-points. 

The inner sub-barrel had fibers oriented parallel to the beam axis, and the outer one had fibers 
inclined by 10° (stereo layer). 

Tracks in the recoil detector are constructed from 3D "space-points" in the SSD and the SFT. 
The ID coordinates from the two sides of a SSD sensor or from adjacent parallel and stereo sub- 
barrels of a SFT barrel are combined to form 2D coordinates. For the SSD, this is accomplished 
taking cnergy-dcposition-correlated combinations of ID coordinates, while for the SFT all geo- 
metrically possible combinations are used. The average (sum) of energy deposits for the two ID 
coordinates is associated with each 2D coordinate for the SSD (SFT). Space-points are constructed 
from 2D coordinates using detector positions. Space-point quadruples are obtained as all possible 
combinations of four space-points, one in each of the two SSD sensors and the two SFT barrels, 
sec figure 3. A "geometrical fit" , which only takes into account coordinate information from the 
space-points, is performed for each such track candidate. This fit uses a helical hypothesis includ- 
ing the lcpton-beam axis, and the track candidate is accepted if the \ 2 value is less than 20. This 
generous value for four degrees of freedom provides high efficiency while removing false tracks that 
arise mainly from space-point quadruples. All possible triples or pairs of space-points are similarly 
treated, whereby space-points belonging to already accepted four-space-point tracks are no longer 
considered. The average beam-axis location was determined by fitting space-point quadruples, with 
frequent corrections for beam movement measured independently by beam-position monitors. 

The momentum of each track is refitted including energy deposition in the SSD under the 
assumption that the particle is a proton, taking into account multiple scattering and energy losses 
in active and passive material. If the resulting value of x 2 exceeds 100 (the optimal value chosen 
after detailed Monte Carlo studies) the refit is discarded under the assumption that the particle 
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Figure 4. Top left: momentum resolution versus momentum for proton reconstruction by using only the 
information on the curvature in the magnetic field (circles) and by combining the information on curvature 
with energy deposition in the SSD (squares). Top right: azimuthal-angle resolution versus momentum. 
Bottom: polar-angle resolution versus momentum. 



is not a proton. In this case, the track is discarded if there are no space-points in the SFT, 
otherwise the momentum reconstruction is based on the geometrical fit. Momentum- resolution 
studies were performed based on Monte Carlo data. In figure 4, the resolution of the momentum 
and angle reconstruction is presented for protons. Reasonable momentum resolution for very low 
momenta is achieved by combining the information on the curvature in the magnetic field with 
energy depositions in the SSD. The azimuthal- and polar-angle resolution is about 4mrad and 
lOmrad, respectively, for proton momenta larger than 0.5 GeV, deteriorating for lower momenta 
because of multiple scattering. 

Particle identification in the recoil detector, described in detail in rcf. [29], is not necessary in 
this analysis because a clean selection of recoil protons is already accomplished by kinematic event 
fitting described in the next section. 

4 Event selection 

In this analysis, inclusive ep — > eX events in the Deep-Inelastic Scattering (DIS) regime are se- 
lected by imposing the following kinematic requirements on the identified positron with the largest 
momentum in the event, as calculated from its four- momentum and that of the incident beam 
positron: 1 GeV 2 < Q 2 < lOGeV 2 , W 2 > 9GeV 2 , v < 22GcV, and 0.03 < x B < 0.35, where 
v = (pq)/M p is the energy of the virtual photon in the target-rest frame, and W the invariant mass 
of the 7*p system [30]. This sample of inclusive DIS events is employed for determination of relative 
luminosities of the two beam-helicity states. 



-7- 



Exclusive ep — > ep7 event candidates are selected from the DIS sample by requiring in the 
forward spectrometer the detection of exactly one identified positron in the absence of other charged 
particles and of exactly one signal cluster in the calorimeter not associated with the positron and 
hence signifying a real photon. The cluster is required to represent an energy deposition above 
5 GeV in the calorimeter and above 1 MeV in the preshower detector. Two kinematic constraints 
that were applied in previous Hermes DVCS analyses to reduce background are also applied here 
in order to maintain compatibility and allow direct comparison: i) the polar angle # 7 » 7 between 
the laboratory three-momenta q and q' is limited to be less than 45 mrad, where q and q' arc 
the three- momenta of the virtual and real photon, respectively (sec figure 1); ii) the value of —t is 
limited to be less than 0.7 GeV 2 . Here, —t is calculated without use of either the photon-energy 
measurement or recoil-detector information, under the hypothesis of an exclusive ep — > epj event 
[7]: 

f = -Q 2 -2u(y- v /^ + Q 2 C0S g 7 , 7 ) 

l + Tt( v ~ V" 2 + Q 2 cos6> 7 . 7 ) 

Moreover, the separation in polar angle between the virtual and real photons is required to be larger 
than 5 mrad. This value is determined mainly by the lepton-momentum resolution. 

All exclusive event samples considered in this paper are derived from the data set collected in 
the years 2006/2007 requiring full functionality of the recoil detector. This data set is a subset of 
that selected in ref. [9] without any such requirement. 

A "pure" exclusive event sample is selected by combining information from the recoil detector 
and forward spectrometer in a kinematic event fit. This fit is based on four- momentum conservation 
under the hypothesis of the process ep — > epj. It is performed for every exclusive-event candidate 
by using the three-momenta of the positron and photon measured in the forward spectrometer and 
the proton candidate in the recoil detector. The quantity 

j* iVfit ~meas\2 

*L = E ~ J ( 4 - 2 ) 
is minimized under the four constraints fj from three-momentum conservation and assumed masses: 

/ i (r? t ) rf t ) ...,Tf t ) = ) j = 1,4, (4.3) 

where r™ cas (rf*) are measured (fitted) kinematic parameters of the positron, photon, and the 
proton candidate and <ii are the measurement uncertainties of these parameters. The minimization 
is conveniently performed using penalty terms: 



X 



fVfit „mcas\2 4 \ f . ( r Rt r fi Ml 

u = e (r< 9 +t 'E [/j( y;;; 9)J . ^ 

j= i \ a 3 1 



where crj are the propagated uncertainties of fj and T is a constant number. For sufficiently large 
T (the value of 10 8 is chosen for this analysis), the constraints are automatically satisfied after 
convergence of the minimization procedure. If more than one proton candidate is reconstructed, the 
one is selected that resulted in the smallest x^ in value from the kinematic event fit. The probability 
calculated from ^kin that a particular event satisfied the ep — > epj hypothesis is required to be 
larger than 0.01, a value that is adequate to ensure negligible background contamination. The 
performance of this event selection is studied using an appropriate mixture of simulated signal 
and background events [31, 32] (the simulation is described near the end of this section). Events 
satisfying all other previously mentioned constraints are found to be selected with high efficiency 
(83%) and background contamination less than 0.2%. This performance is clearly superior to that 
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from imposing only individual constraints on, e.g., the difference between the proton-candidate 
azimuthal angle or transverse momentum measured by the recoil detector and the expected value 
of the corresponding variable calculated from the four-momenta of the positron and the real photon 
detected by the forward spectrometer. 

In the analysis of data collected prior to the installation of the recoil detector, and in the 
analysis of the present data set without using recoil-detector information, the selection of exclusive 
ep — > epj events is performed by requiring the square of the missing mass 

M 2 x = {k+p-k' -q>)\ (4.5) 

calculated using the four-momenta of only the lepton and the real photon, to be within an "exclusive 
region" about the squared proton mass, with boundaries defined by the resolution of the forward 
spectrometer: — (1.5GeV) 2 < M\ < (1.7GeV) 2 . Such an event sample includes not only ep 
epj events but also contamination from "associated" production, ep — > eNnj, including resonant 
production ep — > eA + "f. This contamination is regarded as unresolved background that remains 
part of the signal in Hermes DVCS analyses that do not use recoil-detector information. (A 
correction is applied for other background, as described in section 6.) It is estimated using the 
mixture of simulated events to be about 12% on average within the exclusive region, as illustrated 
in figure 6. Such an exclusive event sample selected by imposing constraints only on the lepton and 
photon four-momenta is named "unresolved" in the following. 

In contrast, the analysis of the pure sample, which includes the reconstruction of the recoil pro- 
ton and kinematic event fitting, introduces two entangled modifications - a background-free mea- 
surement and the kinematic restriction imposed by the acceptance of the recoil detector. In order to 
separate these two effects, the results from the pure sample are compared to results from a subset of 
the unresolved sample that is subject to the same kinematic restriction. This "unresolved-reference" 
event sample is selected from the unresolved sample by requiring the missing four-momentum ( "hy- 
pothetical proton" ) to be within the acceptance of the recoil detector. This requirement results 
in a loss of about 24% of the events. One source of the loss is the effect of the gaps between the 
SSD modules. The other main source is loss of recoil protons with p < 125 MeV, i.e., protons that 
have too low a momentum to reach the outer layer of the SSD because they are stopped in cither 
the target cell or in the inner layer of the SSD. This lower momentum threshold corresponds to 
loss of events at low values of —t < 0.016 GcV 2 . Requiring the proton to be in the recoil-detector 
acceptance leads to a small modification of the average values (— t), (Q 2 ), and (xb) in each kine- 
matic bin compared to the values without such a requirement, as shown in table 1. As expected by 
construction of the unresolved-reference sample, the table demonstrates that the average kinematic 
values of this sample are very similar to those of the pure sample, ensuring that the observables for 
exclusive photon production are the same for the two samples. 

Table 2 summarizes the number of collected events for each of the three exclusive samples: 
unresolved, unresolved-reference, pure, and the average values of the lepton-beam polarization Pi. 
The yield of pure events represents about 65% of the unresolved-reference yield. Of the total 
35% loss, according to the Monte Carlo studies, the event selection based on kinematic event fitting 
eliminates from the unresolved-reference sample about 17% of background events. This also removes 
17% of ep — > epj events. The remaining 1-2% is attributed to recoil-detector inefficiencies [22]. 

Figure 5 shows luminosity-normalized distributions in Mj*. (eq. (4.5)) for each of the three 
exclusive samples. The figure also presents a comparison of experimental data to a mixture of 
simulated data samples. Bethe-Heitler events are simulated using the Mo-Tsai formalism [33], by 
an event generator based on ref. [32] and described in detail in ref. [34] . This sample of BH events 
includes events from associated production generated using the parameterization of the form factor 
for the resonance region from ref. [35] . (The DVCS process is not included since an event generator 
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Table 1. Average kinematic values for each bin in which the Fourier amplitudes of the beam-helicity 
asymmetry are extracted, for each of the three exclusive samples. The "overall" bin represents a single 
kinematic bin covering the entire kinematic acceptance of the Hermes apparatus. The quantities —t and 
Q 2 are given in units of GeV 2 . 
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Table 2. Hermes data sets on single-photon production collected with fully commissioned recoil detector 
in the years 2006 and 2007, using an unpolarized hydrogen target and a positron beam with longitudinal 
beam polarization Pi. For each beam helicity separately and for the total data set, the following quantities 
are given: the integrated luminosity with a systematic uncertainty of 16% (not used in this analysis), the 
respective numbers of selected events in the DIS sample used for normalization and in the three exclusive 
samples, and the average values of the beam polarization. This polarization has a total relative uncertainty 
of 1.96% (dominated by the systematic uncertainty). 



for associated production in DVCS is unavailable.) Semi-inclusive events are simulated by using an 
event generator based on LEPTO [31], including the RADGEN [36] package for radiative effects. 
The Monte-Carlo yield exceeds the experimental data by about 20% in the exclusive region, as 
observed in previous studies of Hermes data [37]. The rightmost panel of figure 5 demonstrates 
that the simulation describes the pure sample well enough to validate the negligible estimate of 
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Figure 5. Distributions of the squared missing mass. The histograms are normalized to the number of 
DIS events. Top figure: distributions from experimental data for the three exclusive event samples discussed 
in the text. The requirements applied on the squared missing mass in order to select (only) the unresolved 
and the unresolved-reference samples are indicated as vertical dashed-dotted lines. The exclusive signal 
is expected around the square of the proton mass, indicated as vertical dashed line. Bottom figure, left: 
unresolved; middle: unresolved-reference; right: pure sample. Experimental data, shown as data points 
(uncertainties covered by symbols), are compared to simulated data. In every panel, the contribution from 
BH ep — > ep7 events is indicated as dashed histogram, and the contributions from associated production 
and semi-inclusive background are shown as hatched histograms. The sum of the simulated distributions 
is shown as solid histogram. See text for discussion. 



background in this sample. 

The fractional contributions of the reaction ep — > epr/ and associated processes, determined 
from the aforementioned mixture of simulated signal and background events, are detailed in figure 
6 in each kinematic bin in which the asymmetry amplitudes are extracted. For the pure sample, the 
contribution of the process ep — > epj is found to be close to 100% and the contribution of events 
from associated processes is close to zero in all kinematic bins. In contrast, for the unresolved 
and unresolved-reference sample the contribution of the associated process is on average about 
12% and 14%, respectively, rising with increasing values of — t; i.e., imposing the acceptance of the 
recoil detector on the unresolved sample has little effect on the background fractions. Therefore, 
comparison of the results from the pure and unresolved-reference samples demonstrates the effects 
of elimination of associated background without changing the experimental acceptance. 
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Figure 6. Fractional contributions from the BH process ep — >• ep-y (closed symbols) and the associated 
BH process ep — > eA + j (open symbols), for each of the exclusive samples. The fractional contributions are 
extracted from Monte Carlo simulations and are presented in the same kinematic binning as the asymmetry 
amplitudes in figures 7 and 8. Symbols for the unresolved (unresolved-reference) sample are shifted to the 
left (right) for better visibility. If the points were plotted without such shifts, a difference would only be 
visible in the first —t bin. 



5 Extraction of single-charge beam-helicity asymmetry amplitudes 

Fourier amplitudes of the single-charge beam-helicity asymmetry A\jj(4>] ) of eq. (2.9) are ex- 
tracted from each of the samples described in section 4. The extraction formalism is described 
in more detail in ref. [37]. It is based on a maximum-likelihood technique [38], which provides a 
bin-free fit in the azimuthal angle (f>. Event weights are employed in the fit in order to account for 
luminosity imbalances with respect to beam polarization. 

Based on eq. (2.7), the distribution of the expectation value of the yield for scattering of a 
longitudinally polarized positron beam from an unpolarized hydrogen target is given by 



(N){ct>\e e ,P i )=£{e e ,P i )T 1 (4>)avv((t>)[l + PiA w (4>-,et)}, (5.1) 

where C denotes the integrated luminosity determined by counting inclusive DIS events and r\ 
the detection efficiency. Here it is assumed that the polarization-dependent cross sections depend 
linearly on the kinematic variables over which the yield is integrated. (A systematic uncertainty 
associated with this assumption is discussed in the next section.) The asymmetry -4lu(0;^) is 
expanded in terms of harmonics in <p in order to extract azimuthal asymmetry amplitudes: 

A W (<P; e t ) c A™ 4, sm ( f> + sin(2</>), (5.2) 

where the approximation is due to the truncation of the infinite Fourier series. Note that j4lu ^ s 
related, but not identical to since there is an additional ^-dependence in the lepton propagators 
in eq. (2.10), and there is another sin0 amplitude sP vcs in eq. (2.10). The former statement also 
holds for an d sf ■ 

As a consistency check for extraneous harmonics caused by the lepton propagators in eq. (2.10) 
and as a test of the normalization of the fit, the maximum likelihood fit was repeated including the 
terms A™^ ^ and A'^ 4 '. As expected, these spurious terms were found to be compatible with zero 
within statistical uncertainties and have negligible impact on the resulting asymmetry amplitudes. 
This provides evidence that the experimental acceptance did not suffer instabilities correlated with 
beam hclicity. 
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6 Background corrections and systematic uncertainties 



6.1 Corrections and uncertainty contributions for the unresolved samples 

The asymmetry amplitudes extracted from the unresolved and unresolved-reference samples are 
corrected for the presence of background that involves semi-inclusive or hard-exclusive neutral 
pseudo-scalar meson production (mainly 7T ), where one of the two photons from the meson decay 
escapes the acceptance of the calorimeter, or the two photons are registered as one calorimeter 
cluster, thus faking a single-photon event candidate. In order to correct the measured amplitude 
^4meas for these background processes, the following procedure is applied in every kinematic bin to 
obtain the asymmetry amplitude Ag na i corrected for background: 

, ^mcas /scmi^scmi /cxcl^cxcl / n \ 

^final = T—j. —? , (6.1) 

^ /semi Jexcl 

where / se mi (/exci) is the fraction of semi- inclusive (exclusive) 7r° events in the data sample and A serm 
(A exc i) is the corresponding asymmetry amplitude. These background fractions are determined from 
Monte Carlo simulations. 

For the estimate of the semi-inclusive background fraction, the event generator LEPTO [31] 
is employed, yielding about 2.7% for the unresolved sample and about 3.1% for the unresolved- 
reference sample, with only weak kinematic dependence. The background asymmetry amplitude 
^4scmi is extracted from experimental data. Neutral pions are reconstructed from a sample of events 
where two photons are detected by requiring the invariant mass of the two photons to be close to 
the mass of the neutral pion: 0.1 GeV< M 77 < 0.17 GcV. In addition, the fractional energy carried 
by the neutral pions is required to be large, E^o/v > 0.8, as only these contribute to the exclusive 
region according to simulations [39]. Simulations showed that the extracted 7r° asymmetry does 
not depend on whether only one or both photons are in the acceptance. Asymmetry amplitudes 
are extracted from the resulting two-photon data sample using the same fit function as that used 
to extract the asymmetry amplitudes for the exclusive ep — > epry measurements. 

Hard-exclusive neutral pseudo-scalar meson production was found to be undetectable at Her- 
mes without using the recoil detector [40], and its isolation is found to be difficult even with its 
use. Therefore, the asymmetry is assumed to be equally probable over the range from -1 to 1 and a 
value of 0± is assigned to the background asymmetry amplitude A exc \. For the estimate of the 
fraction of exclusive ir° events, an exclusive event generator is used that is described in greater detail 
in ref. [41]. For both unresolved samples, the contribution from hard-exclusive meson production 
/cxci is about 0.4%. 

The systematic uncertainty for a given corrected amplitude due to the background (bg) con- 
tamination is taken as one half of the correction to the amplitude: SA^ st = ^|^4fi na i — ^4 m cas|- This 
is considered to be sufficient to account for the mismatch of measured and simulated shape (and 
size) of the semi- inclusive background observed in figure 5. The statistical uncertainty on the total 
correction due to the estimated background asymmetry amplitudes A scm i and A exc i is propagated 
as a contribution to the final statistical uncertainty. 

After applying cq. (6.1), the asymmetry amplitude Afi na i extracted from the unresolved or 
unresolved-reference event samples is expected to originate only from ep — > epj and associated 
processes. No correction is applied for the contamination by these associated processes because the 
relevant asymmetry amplitudes are not known, i.e., for these samples the contribution of associated 
production is considered to be part of the signal (see also figure 6). 

6.2 Corrections and uncertainty contributions for the pure sample 

The asymmetry amplitudes extracted from the pure event sample are not corrected for background 
because the estimated contribution to the yield from background processes is negligible: on average 
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about 0.015% from semi-inclusive processes and less than 0.2% from associated processes. 

6.3 Uncertainty contributions common to all samples 

Systematic uncertainties arising from the forward-spectrometer acceptance, smearing, and finite 
bin width have been studied for several previous DVCS analyses (see, e.g., ref. [37]). An "all- 
in-one" estimate of the combination of these uncertainties was designed to account for possible 
discrepancies between an evaluation of model asymmetries at the measured mean kinematic values 
and an evaluation accounting for the experimental limitations listed above using a full simulation of 
the detector. A similar approach is taken here, except that only variant A of ref. [32] is employed, 
which is based on the model of ref. [43] . This variant describes the asymmetry amplitudes extracted 
from the pure sample well, while the other variants employed in ref. [37] do not. Events from 
associated production are not included for this study. Even though the results are plotted versus 
mean values of measured kinematic quantities, the experimental acceptance can influence the results 
because of a non-linear dependence of the observables on those kinematic quantities. The effect 
of the acceptance of the recoil detector on the results from the pure sample is included in this 
study. Because of an improved survey of the apparatus, effects of a misalignment of the forward 
spectrometer are considered to be negligible. Studies have shown that possible misalignments of the 
recoil detector affect only the efficiency of the constraint from the kinematic event fit and not the 
resulting asymmetry amplitudes. The effects of inefficiencies in the trigger scintillator hodoscopes 
of the forward spectrometer and in the recoil detector are found to be negligible. Asymmetry 
results for the pure sample obtained using one of the four quadrants of the recoil detector at a time 
are found to be consistent within statistical uncertainties. No systematic uncertainty was assigned 
for the relative luminosity measurement based on counting inclusive DIS events because possible 
changes in spectrometer-detection efficiencies correlated with beam helicity would tend to cancel in 
the extracted asymmetries, and because the extracted cos(0</)) asymmetry amplitude was found to 
be consistent with zero, as mentioned at the end of section 5. 

6.4 Summary of systematic uncertainties 

The individual contributions to the systematic uncertainty of each of the asymmetry amplitudes 
are added in quadrature to obtain the final systematic uncertainty in every bin in which the am- 
plitudes are extracted. This procedure is applied separately for each of the three exclusive event 
samples. Table 3 gives the results of the background-corrected asymmetry amplitudes extracted in 
a single kinematic bin covering the entire kinematic acceptance of the Hermes apparatus, together 
with their statistical and systematic uncertainties and the individual contributions to the latter. 
The dominant contribution to the systematic uncertainty is the "all-in-one" uncertainty. (This 
contribution is larger than that estimated in refs. [8, 9] because of the use of a different theoretical 
model.) 

There is a separate scale uncertainty of 1.96% arising from the uncertainty in the measurement 
of the beam polarization. This uncertainty is not included in the error bands used to display the 
combined systematic uncertainty. 

7 Results and discussion 

Figure 7 and table 4 show the sin(n^>) amplitudes of the single-charge beam-helicity asymmetry 
extracted from 2006 and 2007 hydrogen data collected with a positron beam and recoil detector. 
The results are displayed in projections versus —t, xb, and Q 2 and also in a single kinematic bin 
covering the entire kinematic acceptance of the Hermes apparatus ( "overall" ) . For the extraction 
and presentation of asymmetry amplitudes, kinematic variables measured by only the forward 
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Table 3. Asymmetry amplitudes extracted from each of the three exclusive event samples in the "overall" 
kinematic bin together with their statistical and systematic uncertainties. Also the individual contributions 
to the latter are given, "bg" indicating the uncertainty arising from the background correction, and "all- 
in-one" the combined uncertainty arising from detector acceptance, smearing, and finite bin width. The 
amplitudes are corrected for background from semi-inclusive and exclusive neutral pions where applicable. 
A separate scale uncertainty arising from the measurement of the beam polarization amounts to 1.96%. 

spectrometer are used. The calculation of xb and Q 2 requires the identification and momentum 
measurement of the scattered lepton. Also for the calculation of —t (see eq. (4.1)), only the measured 
kinematic parameters of the lepton are used. 

7.1 Results for the pure sample 

In figure 7, the results indicated by the circles are extracted from the pure ep — > ep^f event sample 
defined in section 4, i.e., they involve the reconstruction of the recoiling proton and kinematic event 
fitting. The other sets of data points will be discussed below. The overall value of the leading 
sin (j) amplitude is negative and significantly different from zero. Its one-dimensional projections in 
Xb and Q 2 (which are highly correlated) reveal no dependences. There is no clear indication for a 
dependence on — t, although this amplitude is expected to approach zero as — t approaches zero (see 
cq. (2.11)). The overall value of the sub-leading sin(20) amplitude is compatible with zero within 
its total experimental uncertainty. 

7.2 Comparison to results without recoil-proton detection 

In order to demonstrate the effects of contributions by associated processes that were included 
in previous Hermes measurements of the beam-helicity asymmetry, figure 7 also shows the data 
points for the sin(n^)) amplitudes extracted from the unresolved and unresolved-reference samples. 
The comparison is performed in two steps in order to isolate these effects from those of the change 
in the experimental acceptance due to the recoil detector. 

i) The triangles represent data points extracted from the unresolved-reference sample, which 
is obtained by a missing-mass analysis without using recoil-detector information, but requiring 
the hypothetical recoil proton to be in the acceptance of the recoil detector. The comparison 
with the amplitudes extracted from the pure sample demonstrates the change of the measured 
amplitudes arising from only the removal of the events from associated production, because both 
sets of amplitudes are measured within the same acceptance of forward spectrometer combined 
with recoil detector. There is an indication that the overall value of the sin0 amplitude for the 
pure sample, —0.328 ± 0.027 (stat.), is larger in magnitude than that of the unresolved-reference 
sample by 0.054 ± 0.016 (stat.). The calculation of the latter uncertainty is based on the good 
approximation that one sample is a sub-sample of the other. The systematic uncertainties are not 
relevant for this comparison because the dominant contributions are fully correlated. 

ii) The stars represent data points extracted from the unresolved sample, which is obtained by a 
missing-mass analysis without using recoil-detector information. Comparison with the amplitudes 
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Figure 7. Amplitudes of the single-charge beam-helicity asymmetry in deeply virtual Compton scattering 
shown in projections of —t, xb, and Q 2 . The "overall" results shown in the very left panel are extracted in 
a single kinematic bin covering the entire kinematic acceptance. Statistical uncertainties are shown by error 
bars. The bands represent the systematic uncertainties of the amplitudes extracted from the pure sample. A 
separate scale uncertainty arising from the measurement of the beam polarization amounts to 1.96%. Shown 
are amplitudes extracted from a) the pure ep — > epj sample (red circles, shown at their kinematic values) , 
i.e., obtained with recoil- proton reconstruction; b) the unresolved-reference sample (blue triangles, shifted 
to the right for better visibility), i.e., without recoil-proton reconstruction but requiring its four-momentum 
to be in the recoil-detector acceptance; c) the unresolved sample (black stars, shifted to the left for better 
visibility), i.e., without requirements from recoil-detector acceptance and reconstruction. The actually 
reconstructed kinematic values are specified in table 1 for every bin in which the amplitudes are presented. 
The latter two sets of amplitudes are subject to an average contribution of 14% and 12%, respectively, for 
associated processes (see figure 6 for the kinematic dependences). All three sets of amplitudes are extracted 
from the same 2006/2007 positron-beam data set and the results are strongly statistically correlated. 



extracted from the unresolved-reference sample from i) demonstrates the impact of the recoil- 
detector acceptance on the observed amplitudes. These two sets of amplitudes are subject to very 
similar background conditions as discussed near the end of section 4. The overall value of the sin^i 
amplitude extracted from the unresolved-reference sample, — 0.274 ± 0.022 (stat.), is observed to 
be slightly larger in magnitude in comparison with that extracted from the unresolved sample, 
amounting to a difference of 0.024 ± 0.011 (stat.). The calculation of this uncertainty accounts for 
fully correlated data samples. 

As elaborated in section 4, the lower-momentum threshold that arises from imposing the recoil- 
detector acceptance results in a loss of acceptance at low values of —t, which is reflected in the 
larger statistical uncertainty for the amplitude extracted from the unresolved-reference sample in 
the lowest —t bin in figure 7. 

The overall resu lt for the unresolved sample, being representative of previous Hermes 

publications, is —0.250 ± 0.019 (stat.) versus —0.328 ± 0.027 (stat.) for the pure sample. The 
results for the overall sin <f> asymmetry amplitude for the unresolved and pure samples differ by 
0.078 ± 0.019 (stat.), arising from both the acceptance of the recoil detector and the elimination of 
background from associated production. 

The sub-leading sin(20)-amplitude is found to be compatible with zero within total experimental 
uncertainties for all three event samples. 
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Figure 8. Amplitudes of the single-charge beam-helicity asymmetry extracted from the pure ep —¥ ep7 
sample obtained with recoil-proton reconstruction. The amplitudes are presented in projections of —t, 
ib, and Q 2 . The "overall" results shown in the very left panel are extracted in a single kinematic bin 
covering the entire kinematic acceptance. Statistical (systematic) uncertainties are represented by error 
bars (bands). A separate scale uncertainty arising from the measurement of the beam polarization amounts 
to 1.96%. The theoretical models, which are described in the text, are evaluated at the average values of 
the kinematic parameters specified in table 1 (the points are interpolated by straight lines) . The thickness 
of the VGG lines represents the range & va i = l...oo. 



7.3 Comparison with theory 



In figure 8, the asymmetry amplitudes extracted from the pure sample are compared with calcula- 
tions [42] , labeled "VGG Regge" , from the GPD model described in rcf. [43] . Variants of this model 
differ in the t dependence of GPD H. Here, a Regge-inspired ansatz for the t dependence is used. 
The skewness dependence is controlled by the b parameter, where b va \ (6 S ca) is a free parameter for 
the valence (sea) quarks. The result of the model calculation depends only very weakly on the value 
of 6 V ai- For the sea quarks, the skewness-independent variant of the model (6 soa = 00) is consistent 
with the data, while a maximal skewness dependence (6 sca = 1) is disfavored. 

Also shown in the figure are the results from model calculations [44] labeled "KM" . This model 
is a dual representation of GPDs with very weakly entangled skewness and t dependences. The t 
dependence is approximated by a physically motivated Regge dependence. The model is constrained 
by previous measurements at Hermes, Jefferson Lab, and the collider experiments at Hera. The 
fits resulting in the solid curves disregard data from experiments at Hall A at Jefferson Lab (KMlOa), 
while the dashed curves include these data (KMlOb). The KM calculations agree well with the 
extracted leading amplitude. 

The observed difference between the asymmetries extracted from the pure and the unresolved- 
reference samples is qualitatively consistent with that predicted by a model calculation [45] using 
a soft pion theorem based on chiral symmetry and a A(1232)-rcsonancc model using the large N c 
limit to relate the GPDs for A excitation to those for the nucleon ground state. In this comparison, 
it is important to note that the kinematic conditions for this model correspond approximately to 
the third —t bin of the present measurement. 
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kinematic bin 


unresolved 


unresolved-reference 


pure 






overall bin 


-0.250 ± 0.019 ± 0.047 


-0.274 ± 0.022 ± 0.037 


-0.328 ± 0.027 ± 0.045 


0.00 <-t< 0.06 
0.06 < -t < 0.14 
0.14 <-t< 0.30 
0.30 <-t< 0.70 


-0.234 ± 0.029 ± 0.013 
-0.215 ± 0.035 ± 0.010 
-0.305 ± 0.041 ± 0.015 
-0.324 ± 0.063 ± 0.019 


-0.282 ± 0.039 ± 0.008 
-0.210 ± 0.038 ± 0.011 
-0.321 ± 0.044 ± 0.015 
-0.353 ± 0.066 ± 0.020 


-0.342 ± 0.047 ± 0.016 
-0.256 ± 0.045 ± 0.014 
-0.409 ± 0.053 ± 0.015 
-0.336 ± 0.094 ± 0.022 


0.03 < x B < 0.07 
0.07 < x B < 0.10 
0.10 < x B < 0.15 
0.15 < x B < 0.35 


-0.249 ± 0.031 ± 0.065 
-0.254 ± 0.036 ± 0.050 
-0.260 ± 0.040 ± 0.041 
-0.234 ± 0.054 ± 0.030 


-0.304 ± 0.039 ± 0.051 
-0.262 ± 0.041 ± 0.045 
-0.281 ± 0.044 ± 0.038 
-0.200 ± 0.062 ± 0.028 


-0.302 ± 0.047 ± 0.066 
-0.345 ± 0.050 ± 0.053 
-0.391 ± 0.055 ± 0.043 
-0.262 ± 0.075 ± 0.027 


1.00 < Q 2 < 1.50 
1.50 < Q 2 < 2.30 
2.30 < Q 2 < 3.50 
3.50 < Q 2 < 10.00 


-0.239 ± 0.036 ± 0.062 
-0.271 ± 0.034 ± 0.054 
-0.222 ± 0.038 ± 0.045 
-0.268 ± 0.042 ± 0.033 


-0.291 ± 0.045 ± 0.044 
-0.318 ± 0.040 ± 0.044 
-0.223 ± 0.043 ± 0.042 
-0.256 ± 0.047 ± 0.031 


-0.306 ± 0.054 ± 0.057 
-0.346 ± 0.049 ± 0.053 
-0.298 ± 0.053 ± 0.047 
-0.365 ± 0.058 ± 0.032 




<sin(20) 


overall bin 


0.004 ± 0.019 ± 0.004 


0.011 ± 0.021 ± 0.004 


0.014 ± 0.026 ± 0.002 


0.00 <-t< 0.06 
0.06 <-t< 0.14 
0.14 < -t < 0.30 
0.30 <-t< 0.70 


-0.003 ± 0.029 ± 0.006 
0.038 ± 0.034 ± 0.007 
-0.022 ± 0.041 ± 0.009 
-0.023 ± 0.065 ± 0.002 


0.023 ± 0.038 ± 0.008 
0.019 ± 0.037 ± 0.007 
-0.009 ± 0.043 ± 0.009 
-0.018 ± 0.069 ± 0.003 


-0.010 ± 0.048 ± 0.005 
0.060 ± 0.044 ± 0.008 
-0.044 ± 0.052 ± 0.011 
0.073 ± 0.098 ± 0.007 


0.03 < x B < 0.07 
0.07 < x B < 0.10 
0.10 <x B < 0.15 
0.15 <x B < 0.35 


-0.003 ± 0.031 ± 0.005 
0.001 ± 0.036 ± 0.010 
0.052 ± 0.039 ± 0.007 
-0.060 ± 0.052 ± 0.004 


0.007 ± 0.038 ± 0.003 
0.015 ± 0.041 ± 0.012 
0.042 ± 0.043 ± 0.007 
-0.088 ± 0.061 ± 0.004 


-0.003 ± 0.047 ± 0.010 
0.013 ± 0.050 ± 0.012 
0.104 ± 0.054 ± 0.005 
-0.081 ± 0.075 ± 0.000 


1.00 < Q 2 < 1.50 
1.50 < Q 2 < 2.30 
2.30 < Q 2 < 3.50 
3.50 < Q 2 < 10.00 


0.002 ± 0.036 ± 0.010 
-0.005 ± 0.034 ± 0.002 
0.012 ± 0.038 ± 0.005 
0.014 ± 0.042 ± 0.002 


0.015 ± 0.044 ± 0.022 
0.009 ± 0.040 ± 0.002 
0.014 ± 0.042 ± 0.005 
0.000 ± 0.047 ± 0.002 


0.003 ± 0.053 ± 0.024 
0.018 ± 0.049 ± 0.004 
0.009 ± 0.053 ± 0.006 
0.032 ± 0.059 ± 0.000 



Table 4. Table of results for the beam-helicity Fourier amplitudes A ± 8A ata ,t. ± <5^4 S yst. extracted from the 
three exclusive event samples discussed in the text: unresolved sample, unresolved-reference sample, and 
pure sample. The correlations between the sin</> and sin(2</>) amplitudes are, for each of the samples, below 
10%. The average kinematic values for each bin are compiled in table 1. The quantities —t and Q 2 are 
given in units of GeV 2 . A separate scale uncertainty arising from the measurement of the beam polarization 
amounts to 1.96%. 

8 Summary 

Azimuthal amplitudes of the beam-helicity asymmetry measured at Hermes in exclusive produc- 
tion of real photons of longitudinally polarized positrons incident on an unpolarizcd hydrogen target 
are presented. The asymmetry arises from the deeply virtual Compton scattering process and its 
interference with the Bethe-Heitler process and is sensitive primarily to GPD H. Azimuthal ampli- 
tudes of this asymmetry are extracted from three exclusive data samples: from a pure event sample 
selected by using information from lepton, photon, and recoil-proton detection in a kinematic event 
fit, and two unresolved event samples analyzed by means of a missing-mass technique. The pure 
sample consists almost entirely of ep — > epj events selected in a kinematically complete measure- 
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ment, while the unresolved samples are estimated to contain an average 12-14% contribution from 
associated production and also an approximate 3% contribution from semi-inclusive production, the 
latter of which is corrected for. One of the two unresolved samples serves as reference sample for the 
pure sample, disentangling the effects of the removal of associated background and of the reduced 
detector acceptance caused by the employment of the recoil detector. The recoil detector allows the 
elimination of this background to a very high extent. After correcting for the acceptance reduction, 
the removal of associated background from the data sample alone is shown to increase the magni- 
tude of the leading asymmetry amplitude by 0.054 ±0.016 to -0.328 ±0.027 (stat.) ± 0.045 (syst.). 
The leading asymmetry amplitude obtained from the pure sample is well described by recent fits 
to previously published data. Consistency is found with a theoretical model when a variant is 
considered that assumes skewncss independence for sea quarks (the model shows no sensitivity to 
the skewness of valence quarks). 
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